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A two level system is considered which has no static dipole moment, e.g. molecule H2 in its 
ground electronic state. If strong enough external field is applied, it will dynamically distort such a 
system and supply it with time (and field) dependent dipole moment. Although it is impossible to 
do so in the undistorted system which has no coupling to the dipole component of the external field, 
having induced in it a dipole moment, the rotational and vibrational dynamics of such system can 
be manipulated using lasers. In this work, a system is considered in which the external perturbation 
dynamically induces the transition dipole moment between only two distinct levels. The aim of the 
work is to show how the driving pulse of the following form: 

F[t] = F m[t] cos[w[t] t] (1) 

can be analytically designed, that will produce Rabi-like complete population oscillations between 
the two levels. 



PACS numbers: 3.65.Sq 
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I. CALCULATION 
A. Dynamical equation 

Consider a two level system for which in unperturbed state no spectroscopically allowed dipole transition exists 
between its two levels (e.g. molecule i^)- Such system cannot be directly manipulated using the standard methods 
of laser control. 

However, if such a system is exposed to sufficiently intensive external perturbation, time-dependent dipole transition 
moment between its two states can be dynamically induced. This dipole transition moment can then be used to 
manipulate the system's internal (rotational, vibrational, ...) dynamics. 

It can be demonstrated [1| that the dynamics of such two-level system subjected to strong coherent external radiation 
is governed by the following equation: 

d_ fc a [t}\ = dF[t] ^[t] f-i ^$ ^[t] t e* * \ /ca[t]\ 

dt\cf>[t]J dt w af} [t] \_e-* ««,.[*]* -it&Muj a p[t}t) \cp[t]J {> 

In this equation: c a ^[t] are the wave-function expansion coefficients; F[t] is the external field, fiij[t] = is 
the field-induced dipole moment of state or dipole transition moment (/x is the dipole moment operator); uj a p = 
\E a — Ep\ is the absolute value of the instantaneous transition frequency between the two states of the system; and 
s a p = Sign[E a — Ep] is the sign of that transition frequency. 



B. The external held 



The external field is constructed so as to fit the following analytic form: 

F[t] = F m[t] cos [u[t] t] (3) 

Here, Fq is the maximum amplitude of the field achieved throughout the pulse, m[t] is the pulse envelope (0 < 
m[t] < 1) and ui[t] is the time-dependent frequency of the external perturbation. All of these three parameters are 
freely tunable. 

It should be noted that such 'analytic' form needs not be fully analytic - pulse envelope and pulse chirp can be 
purely numerical functions. Also, such seemingly rigid form of the pulse puts no serious restriction on the possible 
variety of actual pulse shape. Indeed, by varying the three pulse parameters, any reasonable pulse shape can be 
constructed. The only important caveat is that pulse must be such that at any time the field varies much more 
rapidly than the envelope. If this is granted, it straightforwardly follows that: 

dF\t] r , r , , , , , 

— » -F m[i\u>[t] sin [w[t]t] 

gj u[t] t _ g— i uj[t] t 

= 1 F m[t] u[t] (4) 



C. Transforming the time variable 



Next, for further analytic calculation convenience, the following transformations and shorthands are introduced: 



dr = F ° m[t] ^fj dt (5) 

A±[r] EE U)[t[T]]±UJ af3 [t[T}] (6) 

Mr] ee 2 u [t ] t sin [u[t{r}] t[r}] ; (i = a, (3) (7) 

a a [t]\ = fca[t[T]]\ ,„s 

ap[t]J ~ {c p [t[r}}) W 
With Eq. (0J included, in terms of these quantities the dynamical equation Eq. J2J becomes: 
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d_ /o q [t] 
dt xaphi 



(9) 

-f a [ T ] s aP [t[T]}(e-' 1 A -[*l 4 - e l A +[*l ^ 



a a t 



s al3 [t[T]]{e' A_[t]t_ e -i .«.„[*] A+wt) _ //3 [ r ] y V^M 

D. Introduction of RWA and Rabi oscillations condition 

Under certain conditions (see e.g. whose validity can be checked retrospectively, rotating wave approximation 
can be introduced to simplify the obtained full dynamical equation of the system, Eq. 0. If RWA is valid, then 
the dynamical significance of the rapidly rotating elements in the dynamical matrix (the ones containing A + [r]) is 
negligible, and these elements can be dropped from all further calculations. Hence, approximate dynamical equation 
is obtained: 



d_ fa a [r}\ _. ( -f a [r] s a0 [t[r}} e~*'^ *\ fa a [ 

dt \a [T]J ~ 1 \s a0 [t[r]] e*"*® A -M * -f [r] ) \ap[ 



(10) 



b a {r)\ _ e -iA(r) fa a (r)\ (V|) 



Now the final unitary transformation is sought 
with: 

ViM o 
v o P2 [ 

where pi[r] and P2[t] are freely adjustable functions, such that the final transformed system vector satisfies: 



A W= o il - (12) 



d fb a [r}\ _,,f0 A fb a [r) 



ctyb [T}J =±i (l 0){b fi [r]) (13) 

As Ea. l|13|l is identical to the standard form of Rabi oscillations equation in a two- level system (for details, see 
e.g. @), the corresponding solution would represent complete population transfer oscillations between the two levels. 
Introducing transformation into i|l(J[l . the following equation is obtained: 

I fb a [r]\ . ( -f a [r] + £pi[t] s a0 [t[r]} e-'^M] a w+^M-^W)) t [r]x a [t a 

dt {b [r]J {s a0 [t[r}} e*(-/»[*M] A-Tr]+<PiM-,»[r]» »M _ //j[t] + _<L p2 [ r] J ^[r]^ 14J 

If this is to be fitted to form 113(1 . the following conditions must be fulfilled: 

^Pi(r) = f a [r], (15) 

|lP s (r) = Mr], (16) 

Sq/3 [t] A_[t] t[r] + ( Pl [r] - p 2 [r}) = 0, (17) 

which can be compactly written as: 

|l(A_[r] t[r]) = -s^M(/aM - //»M) (18) 
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E. Analyticaly optimized frequency chirp 

Integrating Eq. (|18|) and reverting to the original time coordinate t yields the recurrent formal solution for the 
optimized driving frequency u>[i\: 

F Q /"* 

oj[t] = uj a p[t] - s a p[t]— / (/i m [ii] - Vpplti]) m[ti] uj[ti] t 1 sm[w[ii] *i] dt x (19) 

As this is recurrent equation, id does not provide directly the optimized frequency. However, with numerical com- 
putational power nowadays available it should be a fairly simple and quick task to obtain the solution using some a 
rather simple computer iteration scheme. 

Hence, using this solution, fully controlled and complete population oscillations can be induced in the strongly 
perturbed system. 
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